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ABSTRACT 

Dispersion  in  the  wake  of  an  obstacle  is  studied  as  a 
step  towards  the  analysis  of  bio-agent  release  in  urban  envi¬ 
ronments.  Available  experimental  measurements  are  used 
to  compare  in  details  numerical  predictions  obtained  using 
Reynolds-Averaged  Navier-Stokes  (RANS)  techniques  and 
Direct  Numerical  Simulations  (DNS).  In  particular,  turbu¬ 
lent  scalar  fluxes  are  compared  to  assess  the  quality  of  sim¬ 
ple  closure  based  on  eddy  diffusivity  assumptions.  The 
results  obtained  using  different  RANS  closures  show  con¬ 
siderable  spread  in  the  results,  with  satisfactory  agreement 
obtained  by  one  model  which  predicts  very  low  levels  of 
turbulence  in  the  obstacle  wake.  On  the  other  hand,  DNS 
is  fairly  accurate  but  very  sensitive  to  adopted  inflow  con¬ 
ditions.  Finally,  in  terms  of  the  turbulent  fluxes  the  RANS 
models  qualitatively  represent  the  vertical  flux  but  fail  to 
predict  the  stream  wise  components  while  DNS  predictions, 
at  least  for  the  highest  Reynolds  number  considered,  repre¬ 
sent  both  components  with  a  fairly  good  level  of  accuracy. 

1  INTRODUCTION  AND  BACKGROUND 

The  analysis  of  scalar  dispersion  in  turbulent  flows  is  rel¬ 
evant  to  a  broad  range  of  applications,  including  investiga¬ 
tion  of  hazardous  releases  and  reaction  capabilities  to  bio¬ 
terrorism.  In  the  last  fifty  years  experimental  and  numeri¬ 
cal  investigations  of  canonical  flows  (Fackrell  and  Robins, 
1982;  Livescu  et  al.,  2000)  have  enhanced  theoretical  un¬ 
derstanding  and  subsequently  enabled  the  development  of 
simplified,  semi-empirical  models.  Although  these  models 
have  been  successfully  employed  in  simplified  flow  con¬ 
figurations  (Sykes  et  al.,  1984),  prediction  of  scalar  dis¬ 
persion  over  complex,  realistic  geometries  remains  chal¬ 
lenging  especially  because  of  large-scale  unsteady  effects 
which  cannot  be  properly  accounted  for  in  a  simple  phe¬ 
nomenological  framework.  On  the  other  hand,  it  is  gen¬ 
erally  accepted  that  detailed  flow  simulations  (for  example 
Large  Eddy  Simulation,  LES)  provide  accurate  predictions 
of  the  turbulence  dynamics  and  the  scalar  mixing  rates,  thus 
promising  to  enhance  our  ability  to  study  turbulent  disper¬ 
sion  in  realistic  environments.  This  in  turn  could  lead  to 
the  development  of  better  reduced-order-models. 

In  this  paper  we  present  a  summary  of  numerical  experi¬ 


ments  aimed  at  establishing  the  reliability  of  detailed  flow 
simulations  for  scalar  dispersion  predictions  in  complex  ge¬ 
ometries  (Rossi  and  Iaccarino,  2008).  We  initially  present 
a  comprehensive  Reynolds- Averaged  Navier-Stokes  based 
analysis  of  a  reference  complex  flow  (Vinc^ont  et  al.,  2000) 
where  various  one-point  statistical  models  -  namely  k  —  £, 
k  —  co  and  Reynolds  Stress  Transport  (RST)  -  are  compared 
to  the  experimental  measurements.  As  a  second  step,  we 
performed  direct  numerical  simulations  of  the  same  exper¬ 
imental  setup  with  the  objective  of  provide  a  more  detailed 
analysis  of  the  interaction  between  turbulence  structures 
and  scalar  dispersion. 

2  PROBLEM  DESCRIPTION 

The  adopted  test  case  consists  of  the  experimental  setup 
of  Vingont  et  al.  (2000),  where  the  scalar  dispersion  from  a 
line  source  downstream  of  a  two-dimensional  square  obsta¬ 
cle  has  been  investigated.  The  obstacle  was  completely  im¬ 
mersed  in  a  turbulent-boundary  layer  with  an  approximate 
ratio  8 /h^l,  where  8  is  the  boundary  layer  thickness  pre¬ 
viously  measured  without  the  obstacle  in  place  and  h  the 
obstacle  height.  The  experimental  measurements  were  car¬ 
ried  out  in  a  water  channel  and  in  a  wind  tunnel  where  the 
flow  was  characterized  by  a  Reynolds  number  of  700  and 
1500,  based  on  h  and  the  free  stream  velocity  Uoo ,  and  by  a 
Schmidt  number  of  2500  and  10-6  respectively.  The  scalar 
was  injected  in  the  main  flow  through  a  rectangular  slot 
of  0.14 h  width  located  1  h  downstream  of  the  obstacle.  The 
flow  and  the  scalar  field  were  subsequently  measured  at  two 
different  stream  wise  locations:  4  h  and  6h.  Available  mea¬ 
sured  quantities  are  the  following:  mean  velocity  and  mean 
concentration,  stream  wise  and  vertical  turbulence  intensity, 
Reynolds  shear  stress  and  turbulent  scalar  fluxes.  For  the 
details  about  the  measurements  techniques  the  reader  is  re¬ 
ferred  to  the  paper  of  Vingont  et  al.  (2000). 

3  RANS  ANALYSIS 

This  section  summarizes  the  RANS-based  analysis  of 
the  experimental  test  case.  The  flow  governing  equations 
are  solved  using  the  finite-volume  code  FLUENT  and  fur¬ 
ther  employing  three  different  turbulence  models:  the  stan¬ 
dard  k  —  £  model  (Launder  and  Spalding,  1972),  the  stan¬ 
dard  k  —  co  model  of  Wilcox  (1998)  and  a  Reynolds  Stress 
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Transport  model  (Launder  et  al.,  1975).  In  addition,  the 
one-equation  model  of  Wolfstein  (1969)  is  used  to  resolve 
the  viscosity-affected  near  wall  region  of  the  flow  along 
with  the  k  —  e  and  the  RST  closures.  On  the  other  hand,  the 
k  —  CO  model  already  incorporates  modifications  for  low- 
Reynolds-number  effects  and  is  therefore  adopted  in  the 
outer  flow  as  well  as  in  the  near- wall  region.  The  “frozen” 
velocity  field  is  finally  employed  in  the  finite- volume  code 
PS-SOLVER,  desgined  for  the  solution  of  the  scalar  trans¬ 
port  equation  in  complex  geometries  (Rossi,  2006,  2008). 
Therefore,  the  scalar  advection  is  approximated  to  a  pas¬ 
sive  mechanism.  In  this  framework,  turbulent  scalar  fluxes 
are  estimated  by  a  standard  eddy-diffusivity  model  assum¬ 
ing  a  constant  turbulent  Schmidt  number.  The  numerical 
model  employed  in  the  computations  is  based  on  the  finite- 
volume  grid  shown  in  Fig.  1,  which  has  an  overall  reso¬ 
lution  of  394  x  160  in  the  stream  wise  and  span  wise  direc¬ 
tions,  respectively,  while  64  x  64  cells  are  placed  around 
the  obstacle.  Note  that  a  second-order  accurate  discretiza¬ 
tion  is  adopted  for  both  the  flow  and  scalar  governing  equa¬ 
tions,  based  on  upwind-biased  reconstruction  for  advective 
terms  (Barth  and  Jespersen,  1989).  The  vertical  profiles  of 


Figure  1 .  Sketch  of  the  computational  grid  (Note  that  only  half  the 
grid  points  are  shown  for  sake  of  clarity). 


mean  stream  wise  velocity  are  presented  in  Fig.  2.  The  re¬ 
sults  show  that  the  mean  flow  is  fairly  well  predicted  by 
the  RST  model  for  both  the  water  and  the  air  setup.  More¬ 
over,  while  in  both  cases  the  k  —  £  closure  gives  a  strong 
underestimation  of  the  separated  region  downstream  of  the 
obstacle,  the  profiles  predicted  by  the  k  —  co  model  are 
in  better  agreement  with  the  experimental  measurements. 
The  reattachment  lengths  given  by  the  three  different  tur¬ 
bulence  closures  are  summerized  in  Tab.  1 .  As  suggested 
by  mean  velocity  profiles,  the  RST  model  leads  to  the  best 
agreement  with  the  experimental  results.  It  is  worth  not¬ 
ing  that  although  the  air  setup  is  characterized  by  an  higher 
Reynolds  number  compared  to  the  water  setup,  the  same 
value  of  reattachment  length  has  been  reported  in  the  ex¬ 
periments.  The  computed  vertical  profiles  of  mean  scalar 
concentration  are  shown  in  Fig.  3.  The  most  interesting 


(a)  (b) 


Figure  2.  Profiles  of  mean  streamwise  velocity:  (a)  water  setup  at 
x  =  4 h,  (b)  water  setup  at  x  =  6h,  (c)  air  setup  at  x  =  4 h,  (d)  air 

setup  at  v  =  6 h\  ( - )k  —  £,{ - )k  —  CO,  ( - )  RST,  (o)  Vingont 

et  al.  (2000). 


Table  1 .  Estimated  reattachment  lengths. 


setup 

exp 

k  —  £ 

k-co 

RST 

water 

Ih 

4.55  h 

11.34/z 

1.62h 

air 

Ih 

5.2  h 

9.25  h 

5.65  h 

result  is  that  an  accurate  prediction  of  local  mean  veloc¬ 
ity  profiles  does  not  guarantee  a  satisfactory  prediction  of 
the  scalar  concentration;  the  key  component  is  indeed  the 
turbulent  scalar  diffusivity.  This  is  clearly  shown  by  the 
comparison  of  predicted  mean  concentration  profiles  given 
by  the  k  —  £  and  RST  closures  which  have  a  different  shape 
but  are  almost  equivalent  in  terms  of  average  concentration. 
They  result  largely  understimated  in  the  water  setup,  while 
a  closer  agreement  with  the  experiments  is  found  at  x  =  6h. 
Although  the  k  —  co  model  gives  overall  a  better  agreement 
with  measurements,  the  results  for  the  air  setup  suggest  a 
significant  influence  of  the  Reynolds  number  on  prediction 
capability  of  adopted  turbulence  closures. 

The  overall  effect  of  the  estimated  level  of  turbulent  scalar 
diffusivity  is  clearly  shown  by  the  contour  plots  of  mean 
scalar  concetration  downstream  of  the  obstacle  in  Fig.  4.  A 
similar  scenario  to  the  one  described  in  the  work  of  Vingont 
et  al.  (2000)  is  predicted,  where  the  scalar  is  firstly  con¬ 
vened  upstream  and  in  the  vertical  direction  within  the 


(C)  (d) 

Figure  3.  Profiles  of  mean  scalar  concentration:  (a)  water  setup 
at  x  =  4 h,  (b)  water  setup  at  x  =  6h,  (c)  air  setup  at  x  =  4 h,  (d)  air 

setup  at  x  =  6 h\  ( - )k  —  £,( - )k  —  CO,  ( - )  RST,  (o)  Vingont 

et  al.  (2000). 
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(e)  (f) 

Figure  4.  Mean  scalar  concentration  downstream  of  the  obstacle: 
(left-column)  water  setup,  (right-column)  air  setup;  (a,b)  k  —  £,  (c,d) 
k  —  CO,  (e,f)  RST;  contours  range  is  0.001-0.9. 

flow  reversal  region  to  fill  the  separation  bubble  over  the 
top  of  the  body,  before  being  convected  downstream  by  the 
primary  flow.  Therefore,  each  one  of  adopted  turbulence 


models  are  able  to  capture  the  salient  features  of  the  mean 
flow:  the  large  flow  reversal  region  behind  the  body  and  the 
smaller  separation  bubble  above  the  obstacle.  However,  a 
different  size  of  the  scalar  wake  is  highlighted  in  the  con¬ 
tour  plots.  As  expected,  the  k  —  £  and  the  RST  closures 
yield  a  larger  wake,  suggesting  a  stronger  diffusive  flux 
across  the  wake  boundary.  Moreover,  although  the  RST 
model  gives  a  better  representation  of  the  wake  develop¬ 
ment,  both  turbulence  closures  give  rise  to  a  sharp  decrease 
of  the  scalar  concentration  after  the  injection  into  the  main 
flow.  It  is  clear  by  looking  to  the  results  in  Fig.s  4(c)  and 
4(d)  that  this  is  not  the  case  when  the  k  —  co  model  is  em¬ 
ployed,  where  the  scalar  concentration  is  almost  uniform  in 
the  region  of  flow  reversal  close  to  the  body.  Futhermore, 
a  higher  concentration  is  predicted  in  the  shear  layer  above 
the  obstacle,  which  finally  results  in  the  development  of  a 
confined  scalar  wake. 

The  profiles  of  turbulent  scalar  fluxes  at  x  =  4 h  are  finally 
presented  in  Fig.  5.  A  reasonable  agreement  with  the  ex¬ 
perimental  data  is  found  for  the  k  —  CO  model,  which  is  able 
to  partially  match  the  location  and  magnitude  of  local  ex¬ 
trema  of  the  vertical  component,  while  no  one  of  adopted 
turbulence  closures  is  able  to  predict  the  streamwise  flux. 
Since  it  has  been  shown  that  the  k  —  co  model  gives  a  re¬ 
liable  estimation  of  mean  scalar  concentration,  it  might  be 
concluded  that  the  net  contribution  of  the  streamwise  com¬ 
ponent  in  the  scalar  transport  equation  is  negligeable.  It 
is  worth  outlining  that  the  experiments  show  a  very  strong 
negative  streamwise  flux  within  the  shear  layer  at  y/h  ~  2. 
This  region  does  not  correspond  to  negative  velocity  (at 
x/h  =  4  the  reverse  flow  is  confined  to  y/h  <  0.75)  and 
therefore  corresponds  to  an  actual  anti- diffusion  mecha¬ 
nism. 

4  DNS  ANALYSIS 

The  DNS -based  analysis  of  the  experimental  setup  are 
carried  out  with  the  objective  of  providing  a  more  de¬ 
tailed  analysis  of  the  interaction  between  turbulent  struc¬ 
tures  and  scalar  dispersion.  The  numerical  solution  of  the 
Navier- Stokes  and  scalar  transport  equations  is  obtained 
using  second-order  accurate  finite- volume  techniques  and 
a  block-hexahedral  grid.  In  this  framework,  the  presence 
of  the  obstacle  is  easily  managed  by  employing  an  unstruc¬ 
tured  grid  connectivity.  However,  the  most  critical  issue  in 
the  analysis  of  spatially-developing  flows  using  direct  and 
large-eddy  simulations  is  represented  by  inflow  and  outflow 
boundary  conditions.  While  it  is  generally  accepted  that 
convective  boundary  conditions  at  the  outlet  of  the  com¬ 
putational  domain  are  suitable  to  the  analysis  of  incom¬ 
pressible  flows,  it  has  been  clearly  established  that  inflow 
conditions  can  strongly  affect  the  downstream  development 
of  the  flow  field  (Le  and  Moin,  1992).  Therefore,  in  the 
present  analysis  two  different  boundary  settings  are  con- 
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Figure  5.  Profiles  of  turbulent  scalar  fluxes:  (a)  streanwise  com¬ 
ponent  of  water  setup,  (b)  vertical  component  of  water  setup  ,  (c) 
streamwise  component  of  air  setup,  (d)  vertical  component  of  air 

setup;  ( - )  k  —  £,  ( - )  k  —  CO,  ( - )  RST,  (o)  Vingont  et  al. 

(2000). 


sidered: 

1.  Constant  velocity  profile/zero-fluctuations  (Uniform 
inflow) 

2.  Boundary  layer  profile/superimposed  fluctuations 
(Turbulent  inflow) 

In  both  cases  the  flow  is  supposed  homogeneous  in  the 
spanwise  direction,  where  a  periodic  condition  is  applied, 
while  a  slip  condition  is  employed  at  the  upper  boundary  of 
the  computational  domain.  The  spatial  extent  of  the  com¬ 
putational  domain  is  35 h  x  I5h  x  8/z  while  the  overall  grid 
resolution  in  wall  units  is  reported  in  Tab.  2.  The  code 


Table  2.  Spatial  extent  of  the  computational  domain  and  grid  res¬ 
olution  in  wall  units. 


setup 

size 

Ax+ 

Ay+ 

Az+ 

water 

229  x  80  x  32 

1.2-38.2 

1.2-30.6 

9.5 

air 

229  x  80  x  32 

2.3-73.3 

2.3-58.6 

18.2 

CDP  (Mahesh  et  al.,  2004)  is  adopted  in  the  analysis  of  the 
uniform  inflow  condition  while  the  turbulent  inflow  simula¬ 
tion  is  performed  with  the  help  of  the  FLUENT  code.  The 


specification  of  the  turbulent  inflow  is  based  on  the  random 
flow  generation  technique  originally  proposed  by  Kraich- 
nan  (1970)  and  modified  later  by  Smirnov  et  al.  (2001). 
Note  that  both  codes  adopt  the  fractional  step  method  for 
the  pressure- velocity  coupling  and  a  low-order  central  in¬ 
terpolation  scheme  (LO-CD)  for  advective  terms  to  guar¬ 
antee  the  conservation  of  kinetic  energy  and  therefore  the 
numerical  stability  (Felten  and  Lund,  2006). 

At  each  time  step,  the  continuity- satisfying  mass  flow  rates 
are  employed  at  cell  faces  to  solve  the  scalar  transport  equa¬ 
tion.  For  the  grid  resolution  adopted  in  this  initial  valida¬ 
tion  of  the  numerical  technique,  the  very  sharp  gradients 
arising  at  the  edge  of  the  scalar  wake  are  such  that  upwind- 
biased  schemes  cannot  be  avoided  to  achieve  numerical  sta¬ 
bility.  Therefore,  in  the  CDP  code  the  QUICK  scheme  of 
Leonard  (1979)  is  adopted  for  the  dicretization  of  advective 
terms.  In  the  case  of  the  PS-SOLVER,  a  second-order  accu¬ 
rate  upwind  scheme  (SOU)  based  on  the  slope-limiter  tech¬ 
nique  originally  introduced  by  Van  Leer  (1974)  and  sub¬ 
sequently  extended  to  unstructured  grids  by  Barth  &  Jes- 
persen  (1989)  is  employed.  A  null-flux  condition  is  finally 
specified  at  the  solid  boundaries  while  a  uniform  scalar  in¬ 
jection  (i.e.  a  constant  profile  without  superimposed  fluc¬ 
tuations)  is  assumed  at  the  slot. 

An  initial  condition  of  a  zero-velocity  field  and  zero- 
concentration  of  the  contaminant  is  employed.  The  solu¬ 
tion  is  then  advanced  in  time  using  a  non-dimensional  time- 
step  size  A =  uTAt/h  of  1.862  x  10-3  and  7.623  x  10-4 
for  the  water  and  the  air  setup,  respectively,  until  an  approx¬ 
imate  statistically  steady-state  is  reached.  At  this  point,  tur¬ 
bulence  statistics  are  computed  over  a  time  interval  of  186 
uTt/h  for  the  case  of  water  and  90  uTt/h  for  the  air  setup. 

4.1  Uniform  inflow 

In  this  section  the  results  obtained  from  the  analysis  of 
the  uniform  inflow  condition  are  presented.  For  sake  of 
completeness  the  computations  have  been  carried  out  using 
both  the  CDP  code  and  the  FLUENT/PS -SOLVER  pack¬ 
age,  while  the  results  previously  obtained  using  the  k  —  co 
model  are  also  reported  to  allow  a  direct  comparison.  Note 
that  the  RANS -based  analysis  has  been  performed  by  spec¬ 
ifying  a  boundary  layer  profile  at  the  inflow. 

The  vertical  profiles  of  mean  streamwise  velocity  are 
shown  in  Fig.  6.  Although  the  separated  regions  are  clearly 
overstimated  in  both  the  numerical  experiments,  in  the  case 
of  the  water  setup  the  shape  of  the  velocity  profile  is  also 
very  different  from  measurements:  at  x  =  4h  the  flow  is 
nearly  stagnant  within  the  recirculating  region,  while  at 
x  =  6 h  the  backflow  velocity  is  almost  constant  up  to  the  lo¬ 
cation  of  the  shear  layer  (y  «  1 .5h).  The  results  obtained  by 
the  CDP  and  the  FLUENT  code  are  almost  identical.  The 
situation  is  considerably  different  in  the  air  setup,  where 
the  backflow  velocity  is  not  consistently  computed  using 


the  two  codes.  This  result  indicates  that  the  different  nu¬ 
merical  techniques  adopted  in  the  two  codes  affect  the  pre¬ 
dictions  of  turbulent  structures.  It  is  also  interesting  to  note 
that  the  velocity  at  the  edge  of  the  profile  is  largely  affected 
by  the  uniform  inflow,  being  strongly  overpredicted  com¬ 
pared  to  both  the  reference  datasets.  The  computed  mean 


Figure  6.  Profiles  of  mean  streamwise  velocity  for  the  uniform  in¬ 
flow  condition:  (a)  water  setup  a\x  =  4h,  (b)  water  setup  a\x  =  6h, 

(c)  air  setup  a\x  =  4 h,  (d)  air  setup  at  x  =  6 h\  ( - )  k  —  CO,  ( - ) 

FLUENT,  ( - )  CDP,  (o)  Vingont  et  al.  (2000). 


scalar  profiles  are  presented  in  Fig.  7.  It  is  clear  from  the 
analysis  of  the  water  setup  in  Fig.s  7(a)  and  7(b)  that  tur¬ 
bulent  dispersion  is  nearly  absent  at  both  streamwise  loca¬ 
tions.  This  is  suggested  by  the  local  peak  of  scalar  con¬ 
centration  within  the  shear-layer  (y  «  1.5 h  —  2 h),  which  is 
very  similar  to  the  laminar  profile  that  would  be  obtained 
in  the  limit  of  Sc,Sct  — >  °o  using  RANS  models.  There¬ 
fore,  in  the  case  of  the  water  setup  the  scalar  dispersion  is 
dominated  by  the  mean  flow  transport.  This  is  also  con¬ 
firmed  by  the  constant  scalar  concentration  within  the  ob¬ 
stacle  height  (y  <  1  h)  where  the  flow  is  nearly  stagnant  (see 
Fig.s  6(a)  and  6(b)).  The  concentration  profiles  also  show 
the  more  diffusive  character  of  the  upwind  discretization 
adopted  in  the  code  PS-SOLVER  compared  to  the  QUICK 
scheme  employed  in  the  CDP  code.  In  the  air  setup,  the 
concentration  is  finally  underpredicted  (see  Fig.s  7(c)  and 
7(d))  and  the  absence  of  high  concentration  in  the  shear- 
layer  indicates  the  presence  of  a  stronger  turbulent  scalar 


flux  in  the  vertical  direction. 
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Figure  7.  Profiles  of  mean  scalar  concentration  for  the  uniform 
inflow  condition:  (a)  water  setup  a\x  =  4 h,  (b)  water  setup  at  x  = 

6 h,  (c)  air  setup  a\x  =  4h,  (d)  air  setup  a\x  =  6h\  ( - )  k  —  CO, 

(-•-)  FLUENT,  ( - )  CDP,  (o)  Vingont  et  al.  (2000). 


4.2  Turbulent  inflow 

The  analysis  of  the  results  obtained  by  adopting  the  uni¬ 
form  inflow  clearly  suggests  a  strong  influence  on  the  de¬ 
velopment  of  the  scalar  wake,  determined  by  the  onset  of 
transition  in  the  shear-layer  developing  from  the  leading 
edge  of  the  square  obstacle.  Therefore,  in  this  section  the 
analysis  is  focused  on  the  effect  of  the  turbulent  inflow  con¬ 
dition  described  at  the  beginning  of  Section  4.  The  com¬ 
putations  are  performed  using  the  FLUENT/PS -SOLVER 
package  and  the  results  will  be  directly  compared  to  those 
obtained  from  the  uniform  inflow  analysis. 

A  first  insight  into  the  effect  of  introducing  randomly  gen¬ 
erated  turbulent  fluctuations  at  the  inflow  is  given  in  Fig.  8. 
In  the  case  of  the  water  setup  the  breakdown  of  the  shear- 
layer  occurs  closer  to  the  obstacle;  the  region  above  the 
scalar  source  is  characterized  by  larger  fluctuations.  On 
the  other  hand  in  the  air  setup  the  opposite  behavior  is  ob¬ 
served,  with  the  shear-layer  becoming  more  regular.  In  the 
water  setup,  the  boundary  layer  approaching  the  obstacle 
changes  the  inception  of  the  shear  layer  and,  as  it  has  been 
shown  by  Zhuang  (1999),  produces  an  increased  stream¬ 
lines  curvature  which  is  destabilizing.  On  the  other  hand 


in  the  air  setup  the  velocity  profile  at  the  inlet  results  in 
a  lower  velocity  at  the  obstacle  height,  and  therefore  in  a 
more  stable  shear  layer.  The  analysis  of  mean  velocity  pro- 
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Figure  8.  Instantaneous  scalar  wake  structure  downstream  of  the 
square  obstacle  in  the  fully  developed  regime:  (a)  water  setup  with 
uniform  inflow,  (b)  water  setup  with  turbulent  inflow,  (c)  air  setup 
with  uniform  inflow,  (d)  air  setup  with  turbulent  inflow. 


files  in  Fig.  9  shows  that  the  turbulent  inflow  condition  has 
a  dominant  effect  on  computed  first-order  statistical  mo¬ 
ments  for  the  water  setup.  The  early  stage  of  transition 
occuring  closer  to  the  obstacle  causes  a  significant  change 
in  the  profiles  shape,  particularly  at  x  =  6h.  Furthermore, 
the  adoption  of  the  boundary  layer  profile  at  the  inlet  of 
the  computational  domain  yields  a  very  large  velocity  re¬ 
duction  at  the  edge  of  the  boundary  layer,  giving  a  closer 
agreement  with  the  experiments.  However,  the  region  of 
flow  reversal  is  still  strongly  overstimated.  In  the  case  of 
the  air  setup  the  present  profiles  confirm  the  early  transi¬ 
tion  of  the  shear-layer  even  in  the  absence  of  fluctuations 
at  the  inlet  section,  the  profiles  shape  being  very  similar 
between  the  uniform  and  turbulent  inflow  conditions.  This 
suggests  that  the  turbulent  transport  is  active  in  both  cases. 
Note  that  as  in  the  case  of  the  water  setup  the  outer  velocity 
in  Fig.s  9(c)  and  9(d)  is  significantly  reduced  by  adopting 
the  boundary  layer  profile  at  the  inflow,  but  the  separated 
region  is  still  overpredicted.  The  analysis  of  the  Reynolds 
stress  uv  in  Fig.  10  shows  that  the  profiles  are  consistent 
with  the  scenario  depicted  for  the  mean  velocity  field.  In 
the  case  of  the  water  setup  the  turbulent  transport  at  v  =  4h 
is  practically  absent  when  the  uniform  inflow  is  employed, 
indicating  that  streamwise  and  vertical  turbulent  fluctua¬ 
tions  are  not  correlated.  If  the  turbulent  inflow  is  activated 
a  significant  contribution  is  present  across  the  shear-layer, 
while  it  is  nearly  zero  up  to  the  obstacle  height.  This  is 
again  in  agreement  with  the  small  change  in  the  mean  ve¬ 
locity  profile  shown  in  Fig.  9(a).  The  flow  conditions  at 
x  =  6h  are  similar;  the  turbulent  transport  is  significantly 
larger  in  the  case  of  the  turbulent  inflow  within  the  obstacle 
size,  resulting  in  the  different  shape  of  the  velocity  profile 
in  Fig.  9(b).  The  Reynolds  stress  profiles  for  the  air  setup 
in  Fig.s  10(c)  and  10(d)  suggest  once  again  the  damping  ef- 
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Figure  9.  Profiles  of  mean  streamwise  velocity:  (a)  water  setup  at 
x  =  4 h,  (b)  water  setup  at  x  =  6h,  (c)  air  setup  at  x  =  4 h,  (d)  air 

setup  at  v  =  6 h\  ( - )  uniform  inflow,  ( - )  turbulent  inflow,  (o) 

Vingont  et  al.  (2000). 


feet  following  from  the  boundary  layer  profile  with  super¬ 
imposed  fluctuations  at  the  inlet.  At  both  streamwise  loca¬ 
tions,  a  significant  reduction  of  the  correlation  between  the 
fluctuations  of  velocity  components  is  found.  The  results 
for  the  mean  scalar  concentration  are  presented  in  Fig.  11. 
The  profiles  shape  for  the  water  setup  are  closely  related 
to  the  mean  velocity  field  shown  in  Fig.s  9(a)  and  9(b).  At 
x  =  4 h  the  effect  of  the  turbulent  inflow  is  not  particularly 
evident  and  the  profile  is  just  shifted  downward  in  the  ver¬ 
tical  direction.  Moreover,  a  sharp  increase  of  the  scalar 
concentration  close  to  y  ~  1.5 h  is  found  using  both  inflow 
conditions.  In  spite  of  that,  the  earlier  onset  of  transition 
in  the  case  of  the  turbulent  inflow  can  be  noted  at  v  =  6 h, 
where  the  mean  scalar  profile  is  in  closer  agreement  with 
the  reference  dataset.  On  the  other  hand,  as  it  can  be  noted 
in  Fig.s  11(c)  and  11(d),  the  overall  effect  of  the  fluctua¬ 
tions  at  the  inlet  is  clear  at  both  the  streamwise  locations  for 
the  air  setup.  The  predicted  profiles  are  consistent  with  the 
reduced  level  of  turbulent  velocity  fluctuations,  resulting  in 
the  higher  scalar  concentration  compared  to  the  uniform  in¬ 
flow  results.  However,  while  the  agreement  with  the  exper¬ 
imental  measurements  is  satisfactory  at  v  =  4 h,  the  profile 
obtained  at  x  =  6h  is  overestimated.  The  analysis  of  turbu¬ 
lent  scalar  fluxes  for  the  water  setup  in  Fig.s  12(a)  and  12(b) 
shows  a  significant  improvement  in  the  computed  profiles 
when  the  turbulent  inflow  is  employed.  Moreover,  while 


Figure  1 0.  Profiles  of  Reynolds  stress  uv\  (a)  water  setup  a\x  = 
4 h,  (b)  water  setup  at  x  =  6h,  (c)  air  setup  at  v  =  4 h,  (d)  air  setup 

atx  =  6 h\  ( - )  uniform  inflow,  ( - )  turbulent  inflow,  (o)  Vingont 

et  al.  (2000). 


Figure  1 1 .  Profiles  of  mean  scalar  concentration:  (a)  water  setup 
at  x  =  4 h,  (b)  water  setup  at  x  =  6 h,  (c)  air  setup  at  x  =  4 h,  (d) 

air  setup  at  v  =  6 h\  ( - )  uniform  inflow,  ( - )  turbulent  inflow,  (o) 

Vingont  et  al.  (2000). 


a  vanishing  streamwise  component  has  been  predicted  by 
RANS -based  models,  both  the  magnitude  and  location  of 
the  negative  peak  are  reasonably  reproduced  by  the  present 
computations.  Nonetheless,  the  weak  positive  peak  is  over¬ 
predicted  and  shifted  upward  in  the  vertical  direction  com¬ 
pared  to  the  experimental  measurements.  The  results  also 
show  that  the  vertical  component  is  far  from  the  reference 
dataset,  while  a  better  agreement  for  turbulent  scalar  fluxes 
is  found  for  the  air  setup  in  Fig.s  12(c)  and  12(d).  This  sce¬ 
nario  is  in  agreement  with  the  results  for  the  mean  scalar 
profile  at  the  same  streamwise  location  (see  Fig.s  9(a)  and 
9(c)),  showing  that  accurate  predictions  of  turbulent  trans¬ 
port  have  a  crucial  role  in  high-fidelity  simulations  of  scalar 
dispersion. 

5  CONCLUSIONS 

A  preliminary  study  of  scalar  dispersion  downstream  of 
a  square  obstacle  has  been  presented  as  a  step  towards  the 
analysis  of  bio-agent  release  in  urban  environments.  The 
comparison  with  available  experimental  measurements  has 
shown  that  the  scalar  concentration  given  by  RANS -based 
models  is  dominated  by  the  predicted  eddy-diffusivity.  Al¬ 
though  the  k  —  co  model  has  been  able  to  give  a  good  agree¬ 
ment  for  the  water  setup,  the  analysis  of  the  air  setup  also 
indicates  that  the  predictive  capability  of  the  model  are 
Reynolds  number  dependent.  On  the  other  hand,  DNS 


yield  greater  reliability  in  the  comparison  with  experimen¬ 
tal  results  but  has  been  found  very  sensitive  to  adopted 
inflow  conditions.  However,  the  results  clearly  show  that 
DNS  are  able  to  properly  evaluate  both  the  streamwise  and 
vertical  components  of  turbulent  scalar  fluxes.  Therefore, 
the  present  work  suggests  that  high-fidelity  simulations  can 
be  a  promising  tool  in  the  analysis  of  scalar  dispersion  in 
complex  geometries. 
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